Emission of correlated photon pairs from superluminal perturbations in dispersive 

media. 
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We develop a perturbative theory that describes a superluminal refractive perturbation propa- 
gating in a dispersive medium and the subsequent excitation of the quantum vacuum zero-point 
fluctuations. We find a process similar to the anomalous Doppler effect: photons are emitted in cor- 
related pairs and mainly within a Cerenkov-like cone, one on the forward and the other in backward 
directions. The number of photon pairs emitted from the perturbation increases strongly with the 
degree of superluminality and under realizable experimental conditions, it can reach up to ~ 1CP 2 
photons per pulse. Moreover, it is in principle possible to engineer the host medium so as to mod- 
ify the effective group refractive index. In the presence of "fast light" media, e.g. a with group 
index smaller than unity, a further ~ 10 x enhancement may be achieved and the photon emission 
spectrum is characterized by two sharp peaks that, in future experiments would clearly identify the 
correlated emission of photon pairs. 
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I. INTRODUCTION 

Since the analysis of Schwinger [l[ a great amount of re- 
search has been devoted to the study of particles produc- 
tion from the vacuum. In the seminal paper Q , Hawking 
predicted that the presence of an horizon around a black 
hole induces the particle production. Actually, the exis- 
tence of the Hawking radiation does not require a gravi- 
tational collapse, but rather the key elements are a quan- 
tum field and an event horizon associated with a curved 
space time metric @~0|- Another mechanism providing 
particle production is the dynamical Casimir effect Q. 
Here, a two dimensional quantum theory of a massless 
scalar field is considered that is influenced by the motion 
of a perfectly reflecting boundary (mirror). The vacuum 
expectation value of the energy-momentum tensor for an 
arbitrary mirror trajectory shows a non vanishing radia- 
tion flux. 

Here we consider a further effect by which particle pro- 
duction is induced by the superluminal motion of a per- 
turbation of the refractive index in a dielectric medium 
with dispersion. Differently from the Hawking effect, the 
production of particles is not induced by the presence 
of an horizon but by the fast motion of the perturba- 
tion. Therefore, this effect can be better recognized in 
the context of the anomalous Doppler effect [9t-[l2ll . Fol- 
lowing the description given by Ginzburg, light emitted 
by a generic moving source will be Doppler shifted ac- 
cording to the formula ui(9) = us'/j\l — v / cn(ui) cos 9\ 
where u/ is the comoving reference frame value of 
the emitted frequency, 9 is the observation angle and 
7 = 1/ y/l — v 2 /c 2 . For v/cn(uj) cos 9 < 1 we have the 
normal Doppler effect but for v/cn(ui) cos 9 > 1, i.e. for 
a superluminal source we have the so-called anomalous 
Doppler effect (see Fig. [IJ . Emission of a photon on be- 
half of the moving source must in general correspond to 



a change in the internal state, e.g. a transition in energy 
from one level to another. In the normal Doppler regime, 
a positive frequency photon will thus be emitted with a 
spontaneous transition from an upper to a lower energy 
state. In much the same way, in the anomalous Doppler 
regime a negative frequency photon (positive frequency 
in the laboratory reference frame) will be spontaneously 
emitted in combination with a transition from a lower 
to a higher state. The energy required for this process 
is provided by the superluminal translational movement. 
The superluminal source will therefore continuously emit 
pairs of photons, one in the anomalous and the other in 
the normal Doppler regions [HI]. Anomalous Doppler 
emission has never actually been observed before due to 
the obvious difficulty in actually realising an experimen- 
tal system in which to observe these effects. 

In this work we extend the same concept proposed to 
generate event horizons in optical media: an intense laser 
pulse focused into the medium will induce a local increase 
in the refractive index and this create a perturbation in 
the dielectric medium (DP) that travels locked with the 
laser pulse itself [!, 0] ■ Laser pulses will typically travel 
at a well-defined group velocity so that this kind of DP is 
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FIG. 1: Schematical representation of normal and anomalous 
Doppler emission from a "source" moving with velocity v. 
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ideal for creating the analogue of an event horizon. How- 
ever, by properly choosing how the laser pulse overlaps 
with the dielectric medium it is possible to generate a 
DP that travels with any arbitrary velocity. For example, 
recent experimental results clearly show that it is possi- 
ble to generate isolated peaks within the laser pulse that 
travel faster than the vacuum speed of light, c 0[l5j|. 
Within the context of these findings we therefore consider 
the problem of how a superluminal DP interacts with the 
zero-point fluctuations of the surrounding quantum vac- 
uum. This same problem was first tackled in Ref. [l3[ 
where it was shown that, in a dispersion-less medium a 
superluminal DP will spontaneously emit photons that 
are excited from the vacuum state. The photons appear 
as correlated pairs and their wavelength decreases from 
the far-infrared to the visible region as the DP velocity 
is increased from ~ c to many times c. Here we extend 
these results to a dispersive medium the results and we 
consider in more detail the nature of the correlated pho- 
tons and possible experimental implementation. We also 
consider how the correlated photon pair emission may 
be enhanced by properly engineering the vacuum states, 
i.e. by using media in which the effective group velocity 
at the emitted frequencies is also superluminal, so-called 
"fast light media" [3. 

II. NONDISPERSIVE CASE 

We start by briefly summarising the main results and 
equations for the non-dispersive case treated in Ref. [l3| . 
An intense light pulse induces the perturbation of the re- 
fractive index through the optical nonlinear Kerr effect. 



This perturbation travels across the medium at the same 
speed of the light pulse. The medium was supposed to 
be nondispersive. The authors find that there is no pair 
production unless the pulse velocity v satisfies the con- 
dition v > ^ , where c is the velocity of light in vacuum 
end no is the uniform and constant background refraction 
index. The model analyzed is based on the perturbativc 
approach introduced by Schiitzhold et al. in [l8j and the 
interaction representation for the electromagnetic field 
in the presence of a dielectric constant e(£, t) depending 
on space and time is considered. Moreover, a uniform 
background and a constant value £& = Uq for the di- 
electric constant are assumed. Thus, the disturbance is 
e(x,t) — £5. Following [l8|, the interaction Hamiltonian 

is Hi = £Tl 2 with £ := | ( e (£,t) ^b) an< ^ ^ e nuim:,cr 
of particles created, labeled by the momentum k and the 
polarization fx, is: 

= E +^)i 2 [! - (% • (i) 

where k = (k,u>r) = (k,ck/no), with k = |fc|, e*£ = k/k, 
the subscripts 1 and 2 label the particles of the emit- 
ted couple, and the particular shape of the perturbation 
is contained in the function £. Considering a Gaussian 
behavior for the refraction index n 2 — c(x, t, y, z): 

n 2 = nl + 27i r 1 e-^- v ^ + y 2+z ^ 2 ' 72 (2) 

the expression for the number of created particles 
becomes [11]: 



Nt = 2V 
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[1 -(% 2 -%^ n 



(3) 



The delta function gives a constraint on the possible val- 
ues of the momenta for the emitted particles. 



III. DISPERSIVE CASE 

A. Perturbation 

When considering actual experiments, the previous 
model should only be considered as an approximation of 
a real medium and one has to take into account the pres- 
ence of the dispersion. In this section we generalize the 
model of [l3| for dispersive media. To introduce disper- 
sion in the model we have to refine the theoretic approach 
leading to the expression for the number of emitted par- 
ticles. Let us write the dispersive electromagnetic action 



as 

S=\j E-eEdx\ (4) 

where sE is 

sE{f,t) = e{r,td t )E{f,t). (5) 

We will not specify the operator e now, but we will as- 
sume for it to be a symmetric operator. The unperturbed 
conjugate momentum is then 

n = eE. (6) 
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If we set 



n w (r) 



e~ lult E(f,t)dt, 
e~ lu]t ii(r,t)dt, 



we can write 



1 



- 1 - 1 

E= -n = — I e 

2n J K e(r,w) 



(7) 

(8) 



(9) 



We introduce a perturbation e(r, idt) i-> e(r, i<9 t ) + fc(r — 
vt,idt). To proceed, we note that now the momentum is 



II = + 



(10) 



To describe the Hamiltonian we have to invert this op- 
erator in order to express the electric field in terms of 
the momentum. This can be done perturbatively. Let us 
introduce an order parameter rj so that 

5e(r — vt, id t ) = r]5e(r — vt, id t ). (11) 

Next we write 

E = E a + r\E x + . . . (12) 
where the dots are higher orders in rj. Then, we have 

Tl = eEa+Ti(5eE +£Ei), (13) 

which gives 

II = eE , (14) 

eE x = -5eE - (15) 

We may therefore quantize the theory perturbatively, by 
using the unperturbed momentum operator II and the 
perturbation 

SE = r)E x = —(SsEo) = — (&-ff). (16) 

e e e 

The action is then 

S = i J ft • Edx 4 

2. e 



n.-(fe-ii)^ 

1 f - 1 - , 4 1 f 1 - 1-4 

- / n • -Udx 4 - - / -n • fe-ndx 4 4 
2./ £ 2./ £ £ 



.(17) 



where the dots arc higher order terms and we have used 
the symmetry of £ _1 (and dropped total time deriva- 
tives) . 

Thus, the Hamiltonian density is 



(18) 



where T-Lq is the unperturbed hamiltonian (with disper- 
sion) and 



2 £ £ 



(19) 



is the first order perturbation, which reduces to the 
hamiltonian used in Refs. [TH, [IH in absence of disper- 



sion. 



The perturbative computation is then carried forward 
starting from this expression. 



IV. THE MODEL 

Consider a perturbation of the refractive index moving 
in the x direction with velocity v, i.e. n(Lu)(uj)+Sn(x — vt) 
and assume that the w dependence of the term Sn is 
negligible. Thus, at the first order in 5n the dielectric 
constant is e(w) 2 + 2n(uj)Sn(x — vt). The perturbation 
term (jT9j) may be taken as an operator that we can write 
formally as 

Hi = -^ESeE = —E (n(id t )5n(x - vt)f E, (20) 

where the apex S implies the We yl p rescription for the 
maximally symmetrised operator [ljl- The perturbative 
approach used here may be justified by evaluating the 
main parameters under investigation. We expect to ob- 
serve a similar behaviour to the non-dispersive case for 
which there is a peak of emission for A max = 3er fl3l ]. 
where a is the radius of the gaussian-shaped perturba- 
tion. It is reasonable to expect that the perturbative 
approach works if the measured wavelength A satisfies 
a/X >> ry. For example, if r] is of the order of 0.01 — 0.001 
the perturbative approach is justified for emitted wave- 
lengths well below 100 nm if a ~ 1 fixn. In this sit- 
uation the operator appearing in equation (|20p can be 
approximated by the more tractable one —^ESeE ~ 

—E (n(idt)Sn(x - vt)) S E ~ -E6n(x - vt)n(id t )E, in 
which the operator n(idt) acts just on the electric field E 
and not on the variation of the refractive index 5n(x— vt). 
This perturbation has the same form as the interaction 
Hamiltonian Hi considered in [l~8|, cf. equation (59). In 
the present case Sn(x — vt) plays the role of £(f, t) and 

the electric field of n. These considerations allows us to 
apply the perturbative scheme of [IH . The expression of 
the electric field E is (see fl7l. equation (32), coherently 
adapted to the notation of |18j): 



E(f,t) 



kf.1 



1/2 



1 



n(u))n g (u}) 



kfj. 



-ik-r 



where n g (uj) 



d(n(uj)cu) 



dn(uj) 

i i L — 



(21) 
and k = 



(k,uj) = (fc, cfc/(n(w)), with k = \k\ and eg = k/k. The 
action of the operator n(id t ) on the electric field E is: 



n(id t )E{id t ) J e~ iUkt E u dk 
n(w)E LJ dk, 



(22) 



where E^ are the Fourier modes in the expansion (|21[) 
of the electric field. Thus, the number of particles of 
momentum k\ and polarization [i\ per unit of volume 
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generated by the disturbance Hi turns out to be: where N^ lPil is the number operator. Expanding the 

product of the fields and the number operator in the 
<Nk lfll >=< ip(t ->■ oo)|JVk lA4l |V(* -> 00 ) > bracket with the vacuum state and defining /(fc) := 

= J | d 4 x 2 < m^i)N^Mx 2 )\Q >(23) (#) V2 we 0btain the ex P ression: 



= J d\, J d±x 2 < 0| Yl /( fc i)/( fc 2)/(fc3)/(fc 4 )aj iA(i e £iAii a£ 



eV fit fif fit e*r at eV 

2M2 K2M2 fc/i «M fc 3 /i 3 «3A»3 fe 4 /i 4 «4M4 



fclMl, fc 2 A»2, 
&3M3, fe 4M4 



< 5n(x 1 )on(a; 2 )n(w2)n(w4)e- 4( ' C2+fel)x e l(fe3+fc4)a: |0 >, 
I 



(24) 



where the argument x of <5?i is a four-vector and in our 
particular case it depends just on x and t as a; — vt. We 
indicate with a tilde the four-dimensional Fourier trans- 
form g(k, ^) = J dtdxdydz g(x,t)e lk x ~ i " : o t and using 

the commutation rules for the operators and at we 
obtain: 

AWi = /( fc i) 2 + k 2) 2 f(k 2 ) 2 [n(wi) + n(^ 2 )] 2 

fc2M2 

(25) 

(e r -e r ) 2 . 

The space time dependence of the perturbation drt is of 
the form Sn(t, x, y, z) = Sn(x — vt, y, z). Hence, using the 
variables u = x — vt, w = x + vt, the Fourier transform 
of Sn is 

SnQtj + k 2 ) = —5 (k lx + k 2x - - ( + f 2 ) ) 
v V u W^i) n(u> 2 ) J J 

(26) 



A. Superluminality emission cones 

As for the non dispersive case this last result is very 
meaningful from a physical point of view. The support 
of the delta distribution gives precise conditions for the 
possible states of the emitted particles. First we ana- 
lyze the condition given by the delta function in the non 
dispersive case. Its argument gives the constraint: 

( ftix - —ki) + (k 2x - —k 2 ) = 0. (27) 
\ vn J \ vn J 

This equation implies that whenever k\ x /k\ > c/(vno), 
the second photon must satisfy k 2x /k 2 < c/(vn ), i.e. out 
of every pair of emitted photons, one photon is emitted 



inside the Cerenkov cone #0 = arccos(c/i>no) (anomalous 
Dopplcr region, sec Fig.[T]), and the other is emitted out- 
side the cone, in the normal Dopplcr region. 
In the dispersive case we cannot always identify distinct 
cones of emission for the photons as before. The delta 
function gives the condition: 

f k lx + (k 2x C --k 2 ) = 0. (28) 

V wn(u>x) J \ vn(w 2 ) J 

Similarly to the non-dispersive case, this equation im- 
plies that whenever k\ x jk\ > c/[vn((di)], the sec- 
ond photon must satisfy k 2x /k 2 < c/[vn(uj 2 )]. In 
terms of emission angles for the two photons, we find 
61 < arccos{c/[i;n(a;i)]} and 9 2 > arccos{c/[i>n(u; 2 )]}. 
Therefore, if arccos{c/[wn(wi)]} > arccos{c/[i>n(u; 2 )]}, 




FIG. 2: Schematical representation of correlated photon pair 
emission from a superluminal perturbation moving with ve- 
locity v in the presence of dispersion: the Cerenkov cone is 
split in two and depending on the medium dispersion and 
frequency of emitted photons, regions in which no photons or 
both photons are emitted may appear. For the case of normal 
dispersion (refractive increases with increasing frequency) , the 
figure shows the situation for wi < U2 (a) and cji > U2 (b). 
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the two cones overlap and if arccos{c/[wn(wi)]} < 
arccos{c/[wn(a;2)]} there is a gap between them. In the 
first case there is a region in which both photons can be 
emitted (differently from the non-dispersive case), whilst 
in the second case there is a region in which no photon at 
all can be emitted. The two cases are shown in Figs.[2ja) 
and (b) respectively. The non-dispersive case is obtained 
in the limiting situation in which arccos{c/[i;ro(wi)]} = 
&TCCOs{c/[vn(ui2)] } ■ 



condensed media or ~ 1 ns in gasseous media [20|. This 
situation may therefore be adequately approximated by a 
tanh-like function. We also note that a plasma front will 
locally reduce the refractive index (Sn < 0) as opposed 
to the Kerr effect that will usually increase the refrac- 
tive index (Sn > 0). Yet the result of the calculations 
depends not on the sign of Sn but only on its amplitude, 
velocity and physical dimensions. 
The shape of the perturbation therefore is taken as: 



Gaussian shape of the perturbation 

We now focus attention on the problem of determin- 
ing the actual number of emitted photons, starting from 
the assumption of a Gaussian shape for the refractive 
index perturbation Sn = r^e~ ^ I ^ x ~ v€) +y +z l This 
form for the perturbation is what one may expect when 
a local refractive index variation is induced by a laser 
pulse through the nonlinear Kerr effect. Indeed, laser 
pulses typically have Gaussian-shaped intensity, /, pro- 
file along both longitudinal and transverse coordinates 
so that Sn2l, where 712 is the nonlinear Kerr index, will 
have the same form. A typical value for n, for example in 
fused silica glass, is ~ 0.001 where n 2 ~ 3 • 10~ 16 cm 2 /W 
and we take I ~ 3 • 10 12 W/cm 2 . Inserting this Sn profile 
in (|25p . using — > V J d^k j^s arL d summing over the 
polarisation states u\ and /12 we obtain: 



A k 



0-V"" 2 



(j 1 



v 2 V n 2 (cji)n 2 (wi) 



d 3 k 2 



(29) 



n 2 (uj 2 )n 2 g (uj 2 ) 
' k 2 k 2 + (k 1 -k 2 ) 2 " 



n(wi) + n(w 2 )Y 



h 2 h 2 



S(k lx + k 2x (u)x +cj 2 )) 

v 



Note that in the nondispersive limit, i.e. dn j^ = 0, the 



du) 



expression for the number of emitted particles Nk 1 re- 
duces to equation ©. 



vt 



Sn = 77 tanh ■ 

Ox 

Again we use the variables u — x — vt, w = x + vt, and 
y and z unchanged. Repeating the same computation as 
above we obtain: 



A 1 



ki 



v 2 V 

LJ 2 



2 n 2 (wi)n 2 (wi) 



d 6 k. 



(30) 



n 2 (w 2 )n 2 (uj2 



■ csch 



2 (k ly +k 2 



— {k lx 



k2x) 

n(u)\) + 



k\k 2 + (fcx ■ k 2 ) 2 
h 2 h 2 



n(uj 2 )Y 
1 



S [ h x + k 2x (wi + wa) 



v 



Form these relations for the total emitted photon num- 
bers it is immediately apparent that it is advantageous 
to keep a large a perturbation as possible with each 
of the longitudinal and transverse dimensions a x . y . z 
contributing equally in the multiplicative pre-factor. As 
shown below, it turns out that the tanh profile emits 
roughly half the number of photons of a dimensionally 
similar Gaussian perturbation, indicating that the main 
role is played by the transient switch-on and off of the 
perturbation. We also note the dependence on the group 
refractive index at the emission frequency. As we shall 
discuss below, this allows an additional degree of freedom 
for enhancing or controlling the photon emission process. 



Hyperbolic tangent shape of the perturbation 

We can consider another shape for the perturbation of 
the refractive index. In particular, let us focus on a pro- 
file given by a hyperbolic tangent function. This allows 
us to investigate how the particular geometry of the per- 
turbation affects the maximum of the wavelengths and 
the number of emitted particles. Moreover, a tanh-shape 
may arise if for example the dominant refractive varia- 
tion is related to plasma generated through multiphoton 
or tunnelling ionisation by an intense laser pulse. Intense 
laser pulses may indeed efficiently ionise the medium 
within the first few optical cycles and thus create a very 
steep moving plasma front followed by a nearly constant 
plasma density that will decay (through electron recom- 
bination) over time scales of the order of ~ 100 fs in 



B. Numerical analysis - Gaussian shape 

The non trivial dispersion relation appearing in equa- 
tion (f2U]) makes the dispersive case more intricate to an- 
alyze, even numerically, than the non dispersive one. We 
perform the numerical analysis in fused silica where, for 
completeness we show the full dispersion law given by the 
Scllmcicr relation 1211 : 



n(A) = Re 



^iA 2 
A 2 



/ 2 



^_ 2 A 2 
A 2 



a 3 A 2 



1 1/2 



I 2 

to 



A 2 -I 



(31) 



where 



ai = 0.473115591 
a 2 = 0.631038719 

03 



0.906404498 



h = 0.0129957170 
h = 4.12809220- 10 -3 
l 3 = 98.7685322. 
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k 2 = (k 2x ,k 2y ,0) = (k 2 cos 9 2 ,k 2 sin 2 ,O), then: 
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FIG. 3: (a) log 10 A r A 1 ,A 2 f° r ki x corresponding to the maxi- 



mum, (b) Maxima of log 10 N\ 



for 0i 



iV. 



0, 02 = 

for 0i = 0, 



7T, and 

2 = 7T, 



(7 = ljirn. (c) Maxima of log 10 
and ct = 2/i m. (d) Relation between the wavelengths of the 
emitted couples for f3 increasing from 2, in the right part, to 
20, in the left part of the picture. 



As in the non dispersive case, the delta distribution in 
equation ()29[) gives a constraint that furnishes a relation 
between the momenta of the produced particles. In prin- 
ciple, a component of the momentum, say k 2x , could be 
expressed as a function of the other ones, so that by in- 
tegrating over k 2y and k 2z one could find explicitly the 
function Nj: giving the distribution of the produced par- 
ticles function of the momentum. However, for a 
generic dispersion law, the constraint cannot be analyti- 
cally solved. We therefore study the critical points of the 
integrand of equation (J2H]) that we will call N^fa- This 
is a function of five variables, because the six compo- 
nents of the momenta are not all independent due to the 
constraint given by the delta. The symmetry of the prob- 
lem allows to restrict the analysis to the planes k\ z = 
and k 2z = reducing the number of variables to three. 
This function describes the distribution of emitted cou- 
ples as function of the momenta k\ and k 2 . Its max- 
ima show the momenta of the couple of particles emitted 
with highest probability. Therefore, the effective number 
of emitted particles can be found integrating this dis- 
tribution. We note that the maximum of the function 
Nki,k 2 i s located around 9\ = and 9 2 = n. Indeed, in 
equation (|2"9")) all factors depend on k\ and k 2 by their 

modulus except for the terms \k\ + k 2 \ 2 and [k\ ■ k 2 ) 2 
depending on the relative orientation between the two 
wave vectors. For fixed k\ and k 2 the function 
is maximum for small \ki + k 2 \ 2 and large (k\ ■ k 2 ) 2 . 
Let be k\ = (ki x ,ki y ,0) = (ki cos#i, k\ sin^i, 0) and 



\ki + k 2 \ 2 = k\ + kl 



2k x k 2 cos(6> 2 -6>i) >k\ 
-2kik 2 



(k ■ k') 2 = [ktk 2 cos(6» 2 - 6»i)] 2 < k\k\ 

The two conditions are both satisfied for 6 2 — 6i = it, 
i.e. the two vectors point in opposite directions. The 
condition given by the delta function in equation ([29]) is 
h x + k 2x - f + = that can be rewritten 

as k\ cos 6>i — k 2 cosfi — - , — - /° 2 . = 0, where 9\ is 

the angle between k\ and the x axis and, consequently, 
7r + 9 2 the angle between k 2 and the x axis. As observed 
before, the maximum of equation (|29[) depends only on 
k\, k 2 and k\ ■ k 2 , Thus, at the maxima we have: 



d 

dfx 



N kl ,k 2 (ki, k 2 ,ki- k 2 ) + \g(ki,k 2 , 9i) 



= — A3(fci,fc 2 ,^i) = A(fci - fc 2 )sin6>i = 0, 
d9i 

where A is a Lagrange multiplier. Being k\ and k 2 in 
general not equal and A ^ 0, we therefore have 9\ = 0. 
We also verified numerically that indeed the maximum of 
Nk x M i s a * $1 = ; We visualize this in Figure[3Ja), that 
shows an example (a = 2 and j3 = 5) of the distribution of 
photon numbers in Log scale. The two straight lines cor- 
respond to the condition k\ = k\ x and k 2 = k 2x and the 
curve represents the relation between Ai and X 2 given by 
the delta function in Eq. ([2"9")) . As discussed, these three 
curves intersect at the maximum of implying that 

indeed maximum emission occurs along the propagation 
direction 6\ = 0. 

We perform the numerical analysis for perturbations with 
radius a = 1 /im and a = 2 fj.ni and for increasing 
values of /3. As explained, to give an estimation of 
the number of emitted couples for every single DP we 



(a) cr = l fim 



(b) (7 = 2 fj,m 



max s^zmax 

m m 



P ^lmax ^2max -^Va 1 ,A 2 

m m 



2 2.51e-6 

3 1.93e-6 

4 1.54e-6 

5 1.26e-6 

6 1.08e-6 

7 0.94e-6 

8 0.83e-6 

9 0.75e-6 

10 0.68e-6 
15 0.47e-6 
20 0.36e-6 



4.98e-6 
3.06e-6 
2.17e-6 
1.66e-6 
1.36e-6 
1.15e-6 
0.99e-6 
0.87e-6 
0.78e-6 
0.52e-6 
0.39e-6 



6.13e-7 
6.68e-6 
3.01e-5 
9.35e-5 
2.33e-4 
5.02e-4 
8.73e-4 
1.74-3 
2.91e-3 
2.09e-2 
8.19e-2 



2 3.93e-6 

3 3.50e-6 

4 2.95e-6 

5 2.49e-6 

6 2.13e-6 

7 1.87e-6 

8 1.65e-6 

9 1.48e-6 

10 1.35e-6 
15 0.92e-6 
20 0.70e-6 



7.02e-6 
5.38e-6 
4.11e-6 
3.26e-6 
2.68e-6 
2.69e-6 
1.96e-6 
1.73e-6 
1.54e-6 
1.00e-6 
0.75e-6 



4.14e-8 
2.14e-6 
1.27e-5 
4.28e-5 
l.lle-4 
2.44e-4 
4.80e-4 
8.69e-4 
1.47e-3 
l.lle-2 
4.63e-2 



TABLE I: Maxima of N\ lt \ 2 for increasing /3 and two values 
of a. 
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have to evaluate the integrand of the equation ([29)) , i.e. 
Nk 1 ,k 2 := dN/(dQ,dQ'dkdk'). In the computation the 
value of 6(0) has been approximated with L/(2ir) [see 
Eq. (|26p] . In the following we report all photon numbers 
in function of wavelength rather than wave- vector (thus 
indicated as N\ lt \ 2 ) and integrated over the azimuthal 
angles (pi and tp2 ■ 

Figures [SJb) and (c) show the calculated photon num- 
bers N\ lt \ 2 in Log scale for the value of A2 at which 
emission is maximum and for 77 = 0.001, two different 
values of the perturbation size, a = 1 and 2 /urn respec- 
tively. Figure Eld) shows both wavelengths Ai = Ai max 
and A2 = \2max at which emission is maximum and for 
2 < P < 20. At lowest /3, Ai and A2 have the largest val- 
ues and monotonically decreases and (3 increases. These 
points are also shown for the two values of a consid- 
ered, 1 fj.m (blue dots) and 2 /xm (red dots). We observe 
that the photon emission is peaked at wavelengths Ximax 
that increases with the perturbation diameter. We also 
note that, as in the non dispersive case (l3| . Xi m ax and 
A2max decrease from the mid-infrared region in to the vis- 
ible region as (3, i.e. the perturbation velocity increases. 
However, differently from the non-dispersive case, we un- 
derline that in general Xi max and \2max are different, i.e. 
the photons are correlated but at different wavelengths. 
However, for increasing j3 this difference becomes smaller 

with Xlmax X 2m ax as j3 -> 00. 

Some numerical values of N\ lt \ 2 from these graphs arc 
given in Table [ ^a)| and |(b)| An estimation of the actual 
number of emitted particles can be obtained by integrat- 
ing the function N\ lt \ 2 over A, 9\ and 02- this furnishes, 
for T] = 0.001, cr = l /im, f3 = 20 and a 5 cm propagation 
distance, a number of emitted photons for every DP in 
an angle of 30 degree of ~ 3 • 10 -4 . 



C. Numerical analysis - Hyperbolic tangent shape 

We perform a similar numerical analysis for a hyper- 
bolic tangent shape of the perturbation of the refractive 
index by studying the integrand of equation (|3"0"|) . Adapt- 
ing the previous argument to this case we search its max- 
ima in the direction of the propagation of the DP, i.e. 
for 9 = 0. We assume a perturbation with dimensions 
a x = 1.1 /im and cr y = a z = 1 fim and r\ = 0.001. An es- 
timation similar to the gaussian case, for j3 = 20 and a 5 
cm propagation distance, gives a number of emitted pho- 
tons for every DP in an angle of 30 degree of ~ 1.5 • 10~ 4 . 



D. Fast light 

We now return to a feature outlined earlier on, i.e. the 
dependence of the photon emissivity on the group index 
at the emitted frequencies [see e.g. Eq. (|29j)]. From our 
equations it is clear that by a reduction in the group ve- 
locity at uj\ and/or U2 may greatly influence the actual 
number of emitted photons at both frequencies. Group 



-10 -6 -2 




FIG. 4: (a) log 10 A r A 1 ,A 2 for 6i = 0, 9 2 = tv and for a hy- 
perbolic tangent shape of the perturbation, (b) Correlated 
photon wavelengths for /3 increasing from 2 (top right hand 
point) to 20 (bottom left point). 



velocities exceeding the speed of light have been observed 
in several experiments [22l426l | whereby the medium is ei- 
ther chosen so as to have an absorption resonance close 
to the frequency of interest or is structurally modified, 
e.g. into Bragg grating structures so that the dispersion 
curve is strongly modified (without absorption). This 
property can therefore be exploited to effectively engi- 
neer the (dispersion properties of the) vacuum states, in- 
creasing the number of emitted photons and, moreover, 
providing an additional tool for investigating the corre- 
lation properties of the produced couples. We emulate 
such superluminal group velocities by introducing a sim- 
ple Lorentzian correction on top a background Sellmeier 
relation. We therefore obtain a dispersion law with a nar- 
row peak whose maximum slope is placed in correspon- 
dence to the point of the highest photon emission for a 
medium characterised only by the background dispersion. 
We give an example of such a modified dispersion law for 
the fused silica in Fig. [SJa) and the corresponding group 
index n g (uj)(X(uj)) - A(w)^^ in Fig. GJb). We then 
perform the analysis for a DP of Gaussian shape with 
rj = 0.001, cr = l fim and traveling at j3 = 20 in such a 
medium. In Figure [SJc) we show the distribution of the 
emitted couples in a medium with the background disper- 
sion law and in Figure EJd) the distribution produced by 
the same perturbation in a medium with the Lorentzian 
correction to the dispersion. As a first comment, we ob- 
serve that the number of emitted particles increases by 
a factor ten due to the dependence of Eq. (|2"9"|) on the 
group indices n g (uj\) and n g (uj2). Moreover, the emission 
spectrum is strongly distorted but most importantly, two 
distinct maxima appear. These two peaks are related to 
enhanced emission separately at Ai or at A2. In virtue 
of the nature of the two emitted photons, enhancement 
at one wavelength will necessarily lead to enhancement 
at the correlated photon wavelength. Therefore, observa- 
tion of these two emission peaks in an experiment could 
be considered as evidence of correlated photon emission 
from the superluminal perturbation by measuring only in 
the forward direction. We performed a similar analysis 
for the case of silicon as this is a widely used material 
in optics and waveguide technology where both slow and 
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FIG. 5: Dispersion law (a) and group index (b) for the modi- 
fied dispersion law of the fused silica, (c) log 10 Nx x ,x 2 f° r the 
standard dispersion law of the fused silica, (d) log 10 7VA 1 ,A 2 
for the modified dispersion law of the fused silica. 




FIG. 6: Dispersion law (a) and group index (b) for the mod- 
ified dispersion law of the silicon, (c) log 10 A r A 1 ,A 2 f° r the 
standard dispersion law of the silicon, (d) log 10 N\ 1 ,\ 2 for 
the modified dispersion law. 



fast light may be engineered. We obtain a similar in- 
crease of the emitted particles as shown in Fig. [6] for a 
perturbation of dimension a — 1 /im and velocity /3 = 4. 



E. Conclusions 

Perturbativc analysis of quantum fluctuation excita- 
tion from a travelling refractive index perturbation in- 
dicates that emission of correlated photon pairs occurs 




v=c/cos9 



low-index substrat' 



FIG. 7: Schematical representation of an experimental layout 
for creating a superluminal refractive index perturbation in a 
waveguide. 



only if v > c/n(uj). This effect bears a strong resem- 
blance to the anomalous Dopplcr effect and in this sense 
represent the first detailed analysis of a setting in which 
the effect may actually be observed. The number of emit- 
ted photon pairs is relatively small close to the threshold 
v — c/n(uj) but increases significantly if v is much larger 
(e.g. 10 — 20 x ) larger than c. Our analysis fully accounts 
for material dispersion and, doing so, introduces an in- 
triguing dependence on the group index or group velocity 
of the emitted photons. This dependence and careful en- 
gineering of the host medium may be used to enhance 
photon emission by at least an order of magnitude and 
also provides direct evidence of correlated emission in the 
form of two correlated peaks in the output spectra. 
From our analysis it would seem that experiments are in- 
deed feasible. Refractive index perturbations of the order 
of those used in this work, 77 ~ 0.001, may be obtained 
by focusing laser pulses in fused silica or other nonlinear 
media and higher values have been observed. According 
to our results this would result in a photon pair emission 
rate of ~ 10~ 6 pairs/pulse at (3 = 20 in a 5 mm long 
waveguide, i.e. ~ 1 photon/second if a MHz repetition- 
rate laser is used. These numbers may be increased by an 
order of magnitude or more e.g. by using fast light media 
or higher DP amplitudes. The high velocities used here 
for the perturbation may be obtained experimentally by 
sending an extended (approximately "plane wave" ) laser 
pulse at an angle onto the host medium that could be, 
for example, a waveguide (see Fig. [7]) . The pulse would 
intersect the waveguide with an angle 9 and the Kerr 
effect would create a DP only at the intersection point 
that travels with speed c/ cos 6* (if we approximate the 
refractive index of air ~ 1). Therefore, an incidence an- 
gle of ~ 80 — 85 degrees (i.e. close to normal incidence 
on the waveguide) would guarantee /3 ~ 5 — 20. We note 
that a similar approach has actually been used in pre- 
vious pioneering experiments in which a DC field was 
converted to TeraHertz radiation by a superluminal laser 
pulse excitation of an array of biased capacitors [13] ■ In 
our case we are exciting a superluminal DP in a dielectric 
medium and correlated photon pairs would be generated 
and collected from the two ends of the waveguide. 
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